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Abstract 

We study numerically the effect of on-site Hubbard interaction U between 
two electrons in the quasiperiodic Harper's equation. In the periodic chain 
limit by mapping the problem to that of one electron in two dimensions with a 
diagonal line of impurities of strength U we demonstrate a band of resonance 
two particle pairing states starting from E = U. In the ballistic (metallic) 
regime we show explicitly interaction-assisted extended pairing states and 
multifractal pairing states in the diffusive (critical) regime. We also obtain 
localized pairing states in the gaps and the created subband due to U, whose 
number increases when going to the localized regime, which are responsible for 
reducing the velocity and the diffusion coefficient in the qualitatively similar 
to the non-interacting case ballistic and diffusive dynamics. In the localized 
regime we find propagation enhancement for small U and stronger localization 
for larger U, as in disordered systems. 
PACS numbers: 71.20.+Hk, 71.30.+h 
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I. INTRODUCTION 

Anderson localization! can be also studied in quasiperiodic systems via the Harper's 
equation, which also describes electrons in a square lattice with an added strong magnetic 
fielo&i superconducting networks!, etc. This model presents a very useful alternative to 
the study of one-dimensional (ID) disordered systems, since apart from localization it can 
also display metallic behavior associated with ballistic motion and critical behavior (mo- 
bility edge) with ordinary diffusion&B, somehow mimicing more realistic three-dimensional 
(3-D) disordered systems. In this paper we report results which relate to the problem of 
what happens to the electronic eigenstates and the corresponding quantum dynamics of two 
electrons moving in a quasiperiodic potential in which the interaction between them is taken 
into account. 

The study of two Hubbard interacting particles (THIP) localized by a random potential 
has been pioneered by Shepelyanskyi. This author and produced very interesting 

analytical and numerical work within the Anderson-Hubbard model, which showed weaken- 
ing of Anderson localization, which is known always to be caused by disorder in ID, due 
to the effect of the two particle interaction. This lead to an enhanced propagation effect 
of the interacting electron pair on scales larger than the single-particle localization length. 
The phenomenon of propagation enhancement due to the interaction was also displayed in 
disordered mesoscopic rings threaded by magnetic flux by showing a pairing effect via a few 
/i/2e-periodic, instead of /i/e-periodic, eigenstates0. However, a previous diagonalization 
study for ID disordered systemlll revealed that a few states in the main band show a weak 
propagation enhancement while states with two locally paired electrons are, usually, even 
more localized in the presence of the interaction. Moreover, it was pointed out that for 
large positive U stronger localization occurs when compared to the non-interacting case. In 
a recent work@, it is clearly demonstrated by a tranfer matrix study that no propagation 
enhancement is possible for THIP in an infinite disordered chain at E = 0. 

These, rather conflicting results can be partially understood by the fact that most works, 
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apart from the direct exact diagonalization or dynamical approaches, rely on a mapping 
of the THIP problem to a superimposed banded random matrix ensemble (SBRME). It 
was suggested that if the interaction is expressed in the non-interacting localized basis a 
random band matrix with additional disorder in the matrix diagonal appears (SBRME) 
and the enhancement of the pair localization length £ oc £f, where £1 is the one-particle 
localization length, is easily obtained!. However, the reduction to a SBRME relies on a 
questionable assumption about chaoticity of the non-interacting localized states within £x> 
so that the relevant matrix model could be probably different0. Moreover, since in the one- 
particle localized basis the interaction is relevant only when the two particles are localized 
around positions close to each other, the obtained localization weakening might vanish for 
an infinite chain@. It must be also pointed out that since in most of the previous works only 
the localized case with finite £1 is considered it could be reasonable to expect more dramatic 
pairing effects for extended non-interacting eigenstates which are always overlapping. The 
model considered in this paper allows to study the fate of extended and critical one-electron 
states in the presence of Hubbard interaction. We find a kind of pairing effect for the two 
particle states in the metallic and the critical regimes with the simultaneous appearance 
of localized pairing states in the gaps and the created subband due to the inter actional. 
These localized states reduce the corresponding THIP dynamics although it remains similar 
in nature to the non-interacting case in all three regimes. However, we obtain a weak 
enhancement of propagation in the localized regime when the interaction is switched on, as 
in disordered systems, but even stronger localization is shown to occur for larger U. 

In sections II, III we introduce the Harper-Hubbard model and consider by diagonal- 
ization methods two interacting electrons moving in a tight binding quasiperiodic potential 
of strength A, for various values of the local electron-electron Hubbard interaction U. In 
the absence of the interaction U it is knownH that there are extended states for A < 2, a 
mobility edge for A = A c = 2 and a finite one-electron localization length £1 = l/log(A/2) 
independent of energy for A > 2. In section IV we present our results from the numerical 
diagonalization of the corresponding two interacting electrons Harper-Hubbard equation, by 



showing explicitly extended pairing states in the metallic case A < 2 or multifractal pairing 
states at the mobility edge A = 2. We do not find extended or multifractal pairing states 
for the insulator but a weak propagation enhancement with the simultaneous appearance of 
localized pairing states^. In section V we address the question of the electron localization 
dynamics in the presence of U. The time evolution of a quantum wave packet in the presence 
of interactions shows ballistic motion for A < 2, diffusion for A = 2 and ceases to expand for 
A > 2 as for the non-interacting metallic, critical and localized regimes, respectively^!. In 
the presence of the interaction a decreasing [/-dependent electron velocity and diffusion co- 
efficient due to the appearance of localized states is obtained for the metallic and the critical 
regimes, respectively. In the localized regime larger localization length is found for small U 
although for higher U stronger localization occurs, in agreement with previous results on 
the Anderson-Hubbard mo deli. Finally, in section VI we discuss our results and present 
the conclusions which arise from the present study. 



II. THE HARPER-HUBBARD MODEL 



The Harper-Hubbard tight binding equation for two interacting particles ii 




H = Ys( C n+l,o C n,a + 4, j(T C n+li(7 ) + £ £ A COs(27T</>n)4 >CT C nj(T 
n=l cr n=l a 

n=X 

where <3 n a and c nj(J are the creation and destruction operators for the electron at site n with 
spin cr, A cos(27T0ra) is the potential at site n, with <fi an irrational number usually chosen as 
the golden mean <p = v/ ^~ 1 l and U is the strength of the local Hubbard interaction between 
the two electrons. The Hilbert space can be conveniently divided into one singlet subspace 
with total spin S = and three triplet subspaces with total spin S = 1, S z = 1,0,-1, 
respectively. The three triplet subspaces are energy degenerate and since they permit no 
double occupation the triplet states are not affected by the Hubbard interaction. In a chain 
of N sites the singlet subspace is spaned in the basis of iV(iV + l)/2 spatially symmetric 
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wave functions 

,// w J 75(4,^1,1 + i,T c ni,|)l 0> form^n 2 , 

|^(rii,n 2 ) > s = < v (2) 

[ c Jh,t c JhJ° > for n i = n 2, 
which are antisymmetric with respect to the exchange of the spins and permit double occu- 
pancy. 



III. METHOD OF CALCULATION 

We carried out exact diagonalization of H in the singlet subspace where the Hubbard 
interaction is relevant and found out all the eigenvalues and eigenvectors for finite N sites 
with various A's and Ws. In order to measure the degree of localization for the interact- 
ing electrons we calculate the one-particle spatial extent in the jth two-electron wave 
function vi a0 



N n 



£ (J) = E E \^lnf\J(n l -x 1 Y + {n 2 -x 2 )\ (3) 

ni=l n2=l 

with mean positions 



N ni 

^1,2 = E E |a£UlV,2 ( 4 ) 

of the electrons 1 and 2, where a$ n2 is the normalized coefficient of the wave function in 
the basis of Eq. (0). It must be mentioned that in the way £ is defined it can estimate the 
spatial extend of each electron averaged over the second electron and is related to a quantity 
known as the participation ratioEl. Moreover, £ should correspond to the true localization 
length if the wave functions decay exponentially. Another important quantity used in this 
study is the mean value of the distance between the two electrons in the chain which can be 

calculated for each two-particle wave function via 

N m 

d U) = E E la&jV-nal. (5) 

?li = l U2 = l 

The distance d measures the correlation between the two electrons so that a small d defines 
a pairing two-electron eigenstate, which can be either delocalized in the metallic regime, 
multifractal in the critical regime or localized mostly in the insulating regime A > 2. 



IV. TWO PARTICLE PAIRING STATES 



We diagonalize the Hamiltonian matrix for Fibonacci number chain lengths N, e.g. 
iV = 89 if the rational approximant of = ^ x is ||, so that the potential is periodic 
with period N. In Fig. 1(a) we plot £ and d versus the corresponding eigenvalue E for 
the ID pure A = case with interaction U = 1. The striking characteristic is a band of 
N states, out of the total of 2N 2 two-particle states, which have extremely small distance 
d starting from the energy E = U where d is presisely zero (see Fig. 1(a)). These pairing 
states have one sharp peak at the diagonal line of the plane where n\ = n 2 which implies that 
the particles always stay very close to each other. In Fig. 1(b), (c), (d), (e) we plot some 
characteristic such wave function amplitudes in the plane of the two-electron coordinates 
rii and n 2 where the non-zero amplitudes appear in or very close to the diagonal line. 

Extended pairing states due to the interaction but having a finite width d are also seen 
in Fig. 2 for the ballistic case A = 1 for U — 1. These states are identified from Fig. 
2(a) by plotting in Fig. 2(b), (c) only some states which have small d. The extended 
pairing states in the metallic regime (A < 2) have their number progressively reduced when 
increasing A towards A = 2. Fig. 3(a) accounts for the critical case A = 2 where still a few 
pairing states are seen, such as in Fig. 3(b), with a displayed kind of multifractality along 
the diagonal!. However, apart from extended or multifractal pairing states we also obtain 
another kind of localized pairing states which occur in pairs of almost identical energies and 
similar amplitude distributions. In Fig. 3(c) one such state is shown where a double peaked 
structure is displayed along the diagonal having small d and misleadingly large £ due to our 
definition of £, since such pairing states are strongly localized in two spatial positions along 
the diagonal. The localized pairing states correspond to a physical picture of localization due 
to MotM and they are more frequently encountered in the insulating A > 2 regime. They 
involve tunneling transitions between the two particle localized states spaced at a distance 
proportional to £i apart, having energies that differ by very small amounts. In Fig. 4(a) 
we demonstrate £ and d in the critical regime A = 2 with a higher value of the interaction 
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strength U = 5. In Fig. 4(b), (c) we show two multifractal pairing states. In the plot of 
Fig. 5 we show £ and d for the insulating regime A = 3 with localized pairing states having 
two maxima (Fig. 5(b), (c)) and no extended or multifractal pairing states survive in this 
case. 

We find that localized two-particle states due to the interaction also appear in the 
metallic and the critical regimes. These pairing states have small d and are located either in 
the gaps or in the subband created by the interaction U. They are identified from Fig. 2(a), 
3(a), 4(a), 5(a) and for the critical A = 2 case also in Fig. 6(a), (b), (c) from the integrated 
density of states which is known to be multifractal "devil's staircase" for non-interacting 
electrons. For the THIP four major gaps (plateus) are seen to coexist with smaller gaps 
on all scales. In the created subband for large positive energy the localized pairing states 
due to the effect of the interaction U are clearly seen. It must be emphasised that localized 
pairing states are found for A < 2 only in the presence of finite interaction (U > 0) and the 
results described in this section did not change qualitatively by varying the system size. 

V. TWO PARTICLE DYNAMICS 

The study of the wave packet dynamics provides a global information for the changes 
due to the interaction of all the relevant wave functions. If we put two electrons at the same 
initial site, e.g. the chain center at t = 0, the mean square displacement < (Arr(t)) 2 > for 
each electron at subsequent times t can be calculated from all the singlet eigensolutions of 
Eq. (1) from the variance 

<Ax 2 (t)>=-<n 2 + n 2 2 >= ]T £ 

m=l ri2=l 

where < ... > denotes quantum average and the factor of 2 in the denominator transforms 
< Ax 2 (t) > to correspond to one electron, in order to agree with previous one-electron 
dynamics for U — 0. For the adopted initial condition, in which the two electrons are at 
the same site, only the singlet states with additional on-site energy due to U are allowed. 
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-iEjtAj)* 
a 0,0 



(nl + n 2 2 )/2, 



(6) 
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Alternatively, we have integrated the corresponding two-dimensional equations of motion 
using a Runge-Kutta algorithm, in order to obtain results for much longer chains of N = 
17711. Fig. 7(a) shows the obtained < (Ax(t)) 2 > for A = 1 where is seen that the ballistic 
motion < (Ax(t)) 2 >oc t 2 , remains valid also for finite U but with a reduced velocity. In 
the critical case A = 2 diffusion with < (Ax(t)) 2 >oc t is obtained in Fig. 7(b) reducing 
in magnitude by increasing U, although a tendency for more enhanced propagation is seen 
when U = 5. For the insulator in Fig. 7(c) < (Ax(t)) 2 > shows many oscillations and 
asymptotically reaches larger values for finite {7=1,5 when compared to U = 0, which 
indicates the familiar weakening of localization due to the interactioriHil. However, for 
very large U = 7, 10 the relevant wave functions for the dynamical process lie mostly in the 
subband created above the main band, which corresponds to localized pairing states having 
much shorter localization lengths, and as a result the mean square displacement becomes 
very short (fig. 7(c)) indicating a localizing effect of the interaction^. Therefore, a decrease 
of the degree of localization due to the interaction is demonstrated in the localized regime 
for not too large U, in agreement with the original reported tendency. 

It must be mentioned that our definition of Eq. (6) focuses on the properties of one 
interacting electron and is different from both a\ = \ < (ni+712) 2 > and o 2 _ =< (ni— 712) 2 > 
introducedHJii to examine coherent propagation of two electrons. Our results diplayed in Fig. 
8 for the metallic, critical and the localized cases show a similar behavior of the electron (a 2 ) 
and pair (eri) propagation as well as for the squared pair size (a 2 ). From these results since 
the obtained propagation behavior in the diagonal and its vertical are similar no coherent 
pair propagation can be concluded although a kind of weak pairing can be seen for the metal 
and the insulator where a\ > o 2 _ but no such effect at the critical point where a\ ~ cr 2 . 
Moreover, in the insulating regime for {7 = 1 (Fig. 8(c)) we note that a 2 is below o\. 
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VI. DISCUSSION - CONCLUSIONS 

It can be showni'0 that the interacting electron problem in the periodic A = case can 
be mapped onto an equation for a single electron moving in a two dimensional lattice with 
a line of impurities of energy U along the diagonal. The impurities naturally lead to iV 
resonance states at energies starting from E = U as seen in Fig. 1(a), having amplitude 
only on the impurity sites along the lattice diagonal. In this way extended pairing states 
naturally appear, for example, precisely at E — U a two-particle pairing state can be found 
exactly (Fig. 1(b)) having a constant amplitude on the diagonal and zero elsewhere. It is 
very well known that such resonance extended states can also appear at certain energies in 
non-interacting ID chains with distributed large segments of identical impurities^. 

The quasiperiodicity in addition to the interaction also permits such an exact mapping 
to a single electron equation moving in a two dimensional lattice with a symmetric potential 
A cos(27T0n 1 ) + A cos(27r0n 2 ) at the coordinate n\, 112, which denote the positions of the two 
electrons, in addition to the line of impurities U <5 ni ,n 2 along the diagonal. In analogy with 
the one-dimensional large impurity case where additional small perturbations are known to 
allow the survival of a weak resonant effect at certain energies, remnant of the extended 
states in the absence of perturbation^!, we similarly obtain a kind of pairing states having 
finite but small distance between the two electrons. In this paper we demonstrate by exact 
diagonalization of the THIP Hamiltonian in a finite quasiperiodic Harper's chain such two- 
particle extended or multifractal pairing states due to the interaction for the metal and at the 
critical point. In these regimes we also find localized pairing states due to the interaction in 
the gaps or the created subband, also according to the Mott resonance theory of localization. 
In the localized regime A > 2 we find mostly localized pairing states with short localization 
lengths. 

In disordered systems previous attempts to consider the electron-electron interactions are 
based on perturbation theoriesil'ii or more rigorous solutions for special casesil'il. In this 
paper by a numerical diagonalization study the interaction between two electrons a novel 
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pairing effect for certain delocalized states is found which occurs via extended or multifractal 
pairing states. Localized pairing states are also found due to the interaction in all regimes!!!. 
Our results are confirmed for both repulsive and attractive on-site Hubbard interactions 
and is probably worth noting that in the case of an attractive interaction (U negative) the 
subband created by the interaction lies below the band bottom of the non-interacting case 
so that the ground state is always a localized pairing state for the metal and the insulator. 
The localized pairing states might also affect physical quantities since they can decrease the 
velocity and the diffusion coefficient in the metalic and the critical regimes. Moreover, in 
the corresponding dynamics we demonstrate a tendency for weakening of localization in the 
insulating A > 2 regime for small U but even stronger localization for higher U. 
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FIGURES 



Fig. 1. (a) The electron spatial extend £ (open circles) and the mean distance d (black 
dots) between two electrons as a function of the two-electron wave function energy E for 
the periodic Hubbard chain of A = with interaction strength U — 1. (b),(c),(d),(e) 
Amplitude distributions for extended pairing states with E = U, E = 4.1225, E = 2.5119, 
and E = 4.1173, respectively, as a function of the two electron coordinates n\ and n<i and 
the same parameters as in (a). 

Fig. 2. (a) The £ and d versus the energy E of the two-electron wave functions in the 
ballistic case A = 1 with the interaction strength U — 1. Amplitude distributions (b) for 
extended pairing states E = —0.0013 and (c) for E = —2.5707. 

Fig. 3. (a) The £ and d versus the energy E of the two-electron wave function for the 
diffusive case with A = A c = 2 and the interaction strength U — 1. Amplitude distributions 
(b) for pairing states E = —0.5231 has a multifractal character and (c) for E = 0.8274 is 
localized pairing state in the Mott sense. 

Fig. 4. (a) The £ and d versus the energy E of the two-electron wave function for the 
critical diffusive case A = 2 with U = 5. The states (b) E = 1.3631 and (c) E = 1.3752 
display a multifractal character. 

Fig. 5. (a) The £ and d versus the energy E of the two-electron wave function for the 
insulating localized case A = 3 with U = 1. The (b) E = -0.0646 and (c) E = 0.9083 
correspond to localized pairing states in the Mott sense. 

Fig. 6. (a), (b), (c) The integrated density of states for the critical case A = 2 with 
various values of the interaction strength U . The main gaps correspond to the plateaus. 
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Fig. 7. The mean square displacement < (Ax(t)) 2 > of a wave packet for the approx- 
imant = ^fff with the two electrons initially located at the chain center of long length 
N = 17711 which ensures that the wave does not reach the ends of the chain. The values 
of the interaction strength U are denoted in the figures: (a) log-log plot for the ballistic 
case A = 1, (b) log-log plot for the diffusive case A = A c = 2 and (c) ordinary plot for the 
localized case A = 3. 

Fig. 8. A comparison between the mean square displacements a 2 =< (Ax(t)) 2 >, the pair 
mean square displacement a\ and the pair size a 2 _ with the rest of parameters as in Fig. 7. 
The values of the interaction strength U =1: (a) log-log plot for the ballistic case A = 1, 
(b) log-log plot for the diffusive case A = A c = 2 and (c) ordinary plot for the localized 

case A = 3. 
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